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Abstract. Using CLEO data, we report on the confirmation of the tj^'Sq) ground state of bottomonium in the radiative decay 
Y(3S) —¥ yt]t>- We determine its mass to be M(rjb) = 9391.8 ± 6.6 ± 2.1 MeV, which corresponds to the hyperfine splitting 
AA4/(15) = 68.5 ±6.6 ±2.0 MeV, and the branching fraction ,5§(Y(3S) -> yr] b ) = (7.1 ± 1.8 ± 1.1) x 1(T 4 . These results 
agree with those previously reported by BaBar. 
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INTRODUCTION 

Bottomonium (\bb)) is the preferred system for the study of the qq interaction. Theoretical calculations for bot- 
tomonium are generally more reliable than for other qq systems due to the small relativistic corrections (v/c rs 0.1). 
However, much about the spectrum of bottomonium states remains unknown. Although potential models predict ~ 26 
bottomonium states to lie below the BB threshold (as illustrated in Fig. 1), only 10 of them had been identified since 
the discovery of bottomonium in 1977 [1] until 2008. Most importantly, the ground state r\b{\ 1 Sq), had not been iden- 
tified, despite many searches by CUSB [2] and CLEO [3] at CESR, and ALPEH [4] and DELPHI [5] at LEP, among 
others. 

This situation changed in July 2008, when the BaBar Collaboration announced the discovery of r\ b (\S) in their data 
for 109 million T(3S) [6, 7] in the radiative decay Y(3S) — > yrji,(lS), with a significance of > 10c, as illustrated in 
Fig. 2. We note that the small peak at Ey fa 920 MeV due to the radiative transition to rji,(lS) is not visible until the 
large smooth background is subtracted. BaBar's success is due not simply due to their large data set, more than an 
order of magnitude larger than previously collected at CLEO, but also due a factor 3 greater background suppression 
by the use of a cut on the "thrust angle". The thrust angle is defined as the angle between the signal transition photon 
candidate and the thrust axis of the rest of the event [8]. 

BaBar's discovery of r\ b is exciting and like all such claims, needs to be confirmed by an independent experiment. 
We report on such a measurement at CLEO. 

ANALYSIS OF CLEO DATA 



At CLEO we have an 18 times smaller data set than BaBar (5.9M compared to 109M), so several improvements in 
analysis procedure were needed in order to identify the r\ b signal. The three main improvements were: performing 
a detailed study of the background parameterization; obtaining accurate parameterizations of the photon peak line 
shapes; and using the power of the thrust angle to separate signal and background to perform a joint analysis of the 
full data in three bins of the thrust angle. 

In Fig. 2, we show the inclusive energy spectra for isolated electromagnetic showers in CLEO Y(3S) and T(25) data, 
after contributions from n — > 77 have been rejected. Three features are immediately visible from these spectra. One, 
the spectra are dominated by a larger, smoothly varying background from other bottomonium decays and decays from 
e + e~ — > 7* slighter hadrons. Two, the visible peaks come from the unresolved transitions Xbj(2P, IP) — > 7Y(15'), 
J = 0, 1,2. Lastly, the peaks from e+e' ^srT(IS) and Y(3S,2S) J7] b (\S) (with E 7 (j] b ) w 600 MeV for Y(25) 
and Ey(i~ib) ~ 920 MeV for Y(35)) are expected to be more than an order of magnitude weaker than the %yj(2P, IP) 
peaks, and they reside on the high energy tails of Xbj(2P, IP). 
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FIGURE 1. Spectrum of the bound bottomonium (bb) states. The Ml transitions from the T(n 3 5j) states to the t\i,{\ 1 Sq) are 
indicated by red arrows for the transitions reported in this paper. 
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FIGURE 2. BaBar results for the observation of T]^ [6]. (Left) The observed inclusive photon spctrum. (Right) The background 
subtracted photon spectrum. The peaks, from left to right, are from T(15) ISR, and r\i,. 



Photon Line Shapes: Photon line shapes in crystal calorimeters, like those used by CLEO, are generally parame- 
terized in terms of the Crystal Ball function. This function combines a Gaussian (width a), with a low-energy power 
law tail (parameters a and n). Knowing the shape of the low energy tail is especially important for the %bj peaks, since 
they are important in determing the shape of the background at the energy of the small rj/, signal. These parameters 
must be determined by fitting a background-free photon peak. In our analysis we use two independent methods. In 
one method, we use the observed shape of photons from a given energy from radiative Bhabha events, and in the 
other we use the shape of photons from the excusive decays %b\ (2-P, If) —> yY(lS), Y(1S) — > Both methods give 
consistent results, and once the parameters are determined, they are fixed in subsequent analyses. 

ISR Peak: All the parameters of the ISR photon peak from Y(nS) — > ytsRY(l.S') were fixed. The yield of the ISR 
peak was estimated by extrapolating the observed yield in CLEO data taken at Y(4S), and was then fixed to this value. 

Background: We find that fitting the smooth background is the most crucial component in determing the results 
for the weak rjf, peak. We made several hundred background fits to the data in each of three bins of | cos Ox 1 1 using 
exponential polynomials of various orders (2,3,4), in various energy regions (500-1340 MeV), and with linear and 
logarithmic binning of the data. We found that many of these fit the data acceptably, and took as our final results the 
average of Ey(r\b), SS{X(nS) — > 717/,), and significance for all the good fits (CL> 10%). The r.m.s. variation of these 
fits was then taken as a measure of the systematic uncertainty in the results from this source, ±1 MeV in Ey, ±10% in 
33(r\b), and ±0.4(7 in significance. 




FIGURE 3. CLEO inclusive photon spectra for (left) T(35) decay, and (right) Y(2S), illustrating the broad features described in 
the text. 
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FIGURE 4. (Left) Thrust angle | cos 6j | distribution for Y(3S) data in the expected r\i, signal region, and the expected distribution 
for the r\b signal from MC simulations. (Right) Inclusive photon spectra for the data in the three thrust angle regions, illustrating 
their different signal/background ratios. 



Joint Analysis in Three Bins of |cos07-|: As Fig. 4 (left) shows, the distribution of thrust angle (|cos 9t\) for the 
background-dominated data is strongly peaked in the forward direction, but for the transition photon to the r\h, the 
distribution is expected to be uniform. This happens because the transition photon is uncorrelated with the particles 
produced by the decay of the 77^, whereas the background photons tend to be correlated with the other particles 
produced in the underlying event. We therefore define three different regions of | cos Qj \ ■ region I ( | cos Qj \ =0 — 0.3), 
region II (| cos Qj \ = 0.3 — 0.7), and region III (| cos Qj \ = 0.7 — 1 .0)). Fig. 4 (right) illustrates the different levels of 
signal and background in the three regions. 

Unlike BaBar, we do not reject events in the | cos 0j- 1 > 0.7 region. Instead, we simultaneously fit the spectra for 
each of the three regions, which lets each region contribute to the total result weighted by its individual signal-to- 
background. We call this method the "joint fit". We have analyzed our data by the joint fit method, and also with 
cos 0j I < 0.7, for comparison with BaBar's method. We find that the joint fit method enhances the significance of the 
r\b identification by ~ 1<7. 



JOINT FIT RESULTS 

We show a representative fit for the three bins of thrust angle in Fig. 5. The results of the fit are: N(r]b) = 
23 1 1 ± 546 counts, E y {r\ h ) = 918.6 ± 6.0(stat) MeV, which corresponds to the hyperfine splitting of 68.5 ± 6.6 MeV, 
^(T(35) yrit) = (7.1 ± 1.8(stat)) x 10~ 4 , and significance 4.1a. 

Although we fix V(r]b) = 10 MeV, because the width of rjb is unknown we find that the branching fraction depends 
on the assumed rj h width linearly as J?(Y(3S) — > yr\ b ) x 10 4 = 5.8 +0.13 (r(77/,) in MeV). Our results are dominated 
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FIGURE 5. Background subtracted spectra and representative fit result for T(35) — ¥ yr^lS). The T7/,(1S) peak is clearly visible 
in | cos Qj | region I and II. 
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FIGURE 6. Background subtracted spectra and representative fit result for Y(25) —¥ yilb(lS). The dashed line corresponds to the 
90% CL upper limit given in the text. The r\ b ( 15) is not seen in any of the three regions of | cos Qj \ . 



TABLE 1. Systematic uncertainties in measurement of T(3S) — > yr\b 
photon energy and branching fractions. 



Uncertainty in 



Source 


E Y (MeV) 


,^(T(3S) -» YV„) 


Background (fn, range, binning) 


±1.0 


±10% 


Photon Energy Calibration 


±1.2 




Photon Energy Resolution 


±0.3 


±2% 


CB and Xbj(2P) Parameters 


±0.7 


±8% 


ISR Yield 


±0.4 


±3% 


Photon Reconstruction 




±2% 


N{Y{3S)) 




±2% 


MC Efficiency 




±7% 


Total 


±1.8 


±15% 



TABLE 2. Summary of f] h results from CLEO and BaBar for hyperfine splitting hM h f(\S) bb and 
branching ratio, and comparison with theoretical predictions. 





AM hf (lS) bh , (MeV) 


3§{X{nS) -> yr) b ) x 10 4 


significance 


T(35) -> YVb (CLEO) [9] 


68.5±6.6±2.0 


7.1±1.8±1.1 


4cr 


(BaBar) [6] 


71.4+^ ±2.7 


4.8±0.5±0.6 


> lOo- 


T(2S) -> yilb (CLEO) [9] 




< 8.4 (90% CL) 




(BaBar) [7] 


66.lt 4 .! ±2.0 


3-9±l.ltoi 


3.0o- 


Lattice (UKQCD+HPQCD) [10] 


61 ±14 






(TWQCD) [11] 


70±5 






(Ehrman) [12] 


37±8 






pQCD (various) 


35 - 100 


0.05 - 25 (Y(3S)) 








0.05 - 15 (T(25)) 





by their statistical uncertainties; the systematic uncertainties were determined conservatively and are given in Table 1. 
We further note that although we quote an uncertainty in photon energy calibration of ±1.2 MeV, the energies of the 
ISR peak and the Xbj(2P) centroid agree with their expected energies within ±0.3 MeV. 

The analysis method described for Y(3S) was also used for our Y(2S) data. Because the background near the 
expected Y(2S) — > Ylbi^S) signal (E y « 610 MeV) is approximately 6 times higher than the corresponding Y(3S) (see 
Fig. 3) transition, no evidence for the excitation of r\ b was observed in any |cos07-| bin, as illustrated in Fig. 6. The 
joint analysis led to an upper limit of ^(Y(2S) — > J7]i,) < 8.4 x 10~ 4 , at 90% confidence level. 

We summarize our results and compare them to BaBar's results in Table 2. Both results are in agreement. In 
Table 2 we also list the theoretical predictions for hyperfine splitting and branching ratio. The various potential model 
predictions vary over a wide range. Lattice calculations for the hyperfine splitting are also given [10, 11, 12], which 
generally agree with the experimental results. 
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